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ABSTRACT 

Here we investigate some aspects of stochastic acceleration of ultrarelativistic elec- 
trons by magnetic turbulence. In particular, we discuss the steady-state energy spectra 
of particles undergoing momentum diffusion due to resonant interactions with turbulent 
MHD modes, taking rigorously into account direct energy losses connected with differ- 
ent radiative cooling processes. For the magnetic turbulence we assume a given power 
spectrum of the type W(k) oc k~ q . In contrast to the previous approaches, however, 
we assume a finite range of turbulent wavevectors k, consider a variety of turbulence 
spectral indexes 1 < q < 2, and concentrate on the case of a very inefficient particle 
escape from the acceleration site. We find that for different cooling and injection con- 
ditions, stochastic acceleration processes tend to establish a modified ultrarelativistic 
Maxwellian distribution of radiating particles, with the high-energy exponential cut- 
off shaped by the interplay between cooling and acceleration rates. For example, if the 
timescale for the dominant radiative process scales with the electron momentum as oc p r , 
the resulting electron energy distribution is of the form n e (p) oc p 2 exp [— ^ (p/p eq ) a ], 
where o = 2 — q — r, and p eq is the equilibrium momentum defined by the balance 
between stochastic acceleration and energy losses timescales. We also discuss in more 
detail the synchrotron and inverse-Compton emission spectra produced by such an elec- 
tron energy distribution, taking into account Klein-Nishina effects. We point out that 
the curvature of the high frequency segments of these spectra, even though being pro- 
duced by the same population of electrons, may be substantially different between the 
synchrotron and inverse-Compton components. 
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Introduction 



Stochastic acceleration of ultrarelativistic particles via scatterings by magnetic inhomogeneities 
was the first process discussed in th e context of generation of a power-law energy distribution of 



cosmic rays (jFermi 



1949 



Davislll956l ). Because the characteristic acceleration timescale for a given 
velocity of magnetic inhomogeneities, say Alfven velocity va, is i a cc oc (va/c)~ 2 , the stochastic 
particle acceleration is often referred as a '2nd-order Fermi process'. For commonly occuring non- 
relativistic turbulence, v\ <C c, turbulent acceleration mechanism is often deemed less efficient 
when compared to acceleration by shocks where the rate of momentum change 5p/p ~ w s h/c (hence 
the name lst-order Fermi process). However, here one also needs repeated crossing of the shock 
front by the particles which can come about via scattering by turbulence upstream and downstream 
of the shock. Thus, again the acceleration rate or timescale is determined by the scattering time 
scale. For nonrelativistic turbulence va <C c, relativistic particles p 3> mc 2 , and high-/? or weakly 
magnetized plasma, this time is shorter than the stochastic acceleration time, which may not 
be the case in many astrophisical plasmas. We note that in a relativistic regime, for example, 
lst-o rder Fermi process encoun t ers several difficult i es in accelerating p articles to high energies 



(e.g.. iNiemiec &: Ostrowskill2006l : iNiemiec et al.ll2006l : iLemoine et al.ll2006l ). while at the same time 



stochastic particle energization may play a major role, since velocities of the turbulent modes may be 
high, va ;$ c. And indeed, 2nd-order Fermi processes were being discussed in the context o f different 



astrophysic al sources of high energy radiation and particles, such as accretion d iscs (e.g.. iLiu et al 



2004 



blazars (e.g. iKatarzyhski et alj|20 06b ; 



Stawarz fc Qstrowski 



200 



2; 



2006), clusters of galaxi es (e.g.. iPetrosianl l 2001; Brunctt i fc Lazarianll2007|) . gamma-rav burst s 



(e.g. , IStern Poutanenll2004r). solar flares (e.g..|Petrosian fc Donaghvlll999l ; IPetrosian Liull2004l ) 



Giebels et al 



Stawarz et al 



20071 ) . or extragalactic large-scale jets (e.g. 



2004 ) . We note, that although turbulent acceleration 



is often a process of choice in modeling high energy emission in different objects, and in fact 
there may be some other yet much less understood mechanisms responsible for generation of such 
(like magnetic reconnection) , evidences for the distributed (or in situ) acceleration process taking 
place in several astrop hysical systems are strong (see, e.g. Jjester et alJl200ll ; iKataoka et alJl200a ; 
Hardcastle et al.1 120071 . in the context of extragalactic jets). 



It was pointed out by ISchlickeiseii (|1984l . Il985l ). that continuous (stochastic) acceleration of 
high energy electrons undergoing radiative energy losses tends to establish their ultrarelativistic 
Maxwellian energy distribution, as long as particle escape from the acceleration site is inefficient. 
This analysis concerned a particular case of acceleration timescale independent on the electrons' 
energy, and the dominant synchrotron- type energy losses. Interestingly, very flat (inverted) elec- 
tron spectra of the ultrarelativistic Maxwellian-type — often approximated as a monoenergetic 
electron distribution — were discussed in the context of flat-spectrum radio emission observed 



2001 



from Sgr A* and several active galactic nuclei (see, e.g.. iBeckert Duschllll997l : iBirk et al 
and references therein). More recently, it was proposed that such 'non-standard' electron spec- 
tra can account for striking high-en e rgy X-ray emission of large-scale jets observed by Chandra 
satellite (jStawarz Qstrowskil 120021 ; IStawarz et al.ll2004l ). or correlated X-ray and 7-ray (TeV) 
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emission from several BL Lac objects detected by the modern ground-base Cherenkov Telescopes 
(jKatarzyhski et al.l l2006al ; iGiebels et al.l 120071 ). In addition, it was shown that narrow electron 
spectra, e.g. Maxwellia n distribution, can ex plain properties of extragalactic high brightness tem- 
perature radio sources (jTsang fc Kirk!l2007al lbl). alleviati ng the difficulties ass ociated with the an- 
ticipated by not observed inverse-Compton catastrophe (jOstorero et al.ll2006l ). 

Motivated by these most recent observational and theoretical results, in this paper we in- 
vestigate further some aspects of stochastic acceleration of ultrarelativistic electrons by magnetic 
turbulence. In particular, we discuss steady-state energy spectra of particles undergoing momentum 
diffusion due to resonant interactions with turbulent MHD modes, taking rigorously into account 
direct in situ energy losses connected with different radiative cooling processes. As described in the 
next section § 2, we use the quasilinear approximation for the wave-particle interactions, assuming 
a given power spectrum W(k) oc k~ q for magnetic turbulence within some finite range of turbulent 
wavevector k\ < k < k%, and consider turbulence spectral indexes in the range 1 < q < 2. In 
section § 3 we provide steady-state solutions to the momentum diffusion equation corresponding to 
the case of no particle escape but different cooling and injection conditions. In section § 4 some par- 
ticular solutions are given corresponding to the case of a finite particle escape from the acceleration 
site. In section § 5 we discuss in more details synchrotron and inverse-Compton emission spectra 
of stochastically accelerated electrons, taking into account Klein-Nishina effects. Final discussion 
and conclusions are presented in the last section § 6 of the paper. 



2. General Description 

Let us denote the phase space density of ultrarelativistc particles by f(x,p,t), such that the 
total number of particles is M(t) = f d 3 x J d?p f(x,p,i). Here the position coordinate x and 
the momentum coordinate p are not the position and the momentum of some particular particle, 
but are fixed to the chosen coordinate space, and therefore are independent variables. In the 
case of collisionless plasma, the function f(x,p,t) satisfies the relativistic Vlasov equation with the 
acceleration term being determined by the Lorentz force due to the average plasma electromagnetic 
field acting on particles. This averaged field can be found, in principle, through the Maxwell 
equations, and such an approach would lead to the exact description of the considered system. 
However, due to strongly non-linear character of the resulting equations (and therefore substantial 
complexity of the problem), in most cases an approximate description is of interest. In the 'test 
particle approach', for example, one assumes configuration of electromagnetic field and solves the 
particle kinetic equation to determine particle spectrum. Further simplification can be achieved if 
one assumes presence of only a small-amplitude turbulence (5E, 5B) in addition to the large-scale 
magnetic fielcf] Bq 3> SB, such that the total plasma fields are B = B$ + 5B and E = SE. 



x Due to the expected high conductivity of the plasma one can neglect the large-scale scale electric field, Eo — . 
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In order to find the evolution of the particle distribution function in the phase space under the 
influence of such fluctuating electromagnetic field, it is convenient to consider an ensemble of the 
distribution functions (all equal at some initial time), such that the appropriate ensemble-averaging 
gives (SB) = (5E) = and f(x,p,t) = (f(x,p,t)) + 5f(x,p,t). It can be then shown via the 'quasi- 
linear approximation' of the Vlasov equation th at the ensemble-average of the distrib ution function 



(f(x,p,t)} satisfies the Fokker-Planck equation (jHall &: Sturrock 



1967 



Melrosd ll968H If, in addi- 



tion, the particle distribution function is only slowly varying in space ('diffusion approximation'), 
and the scattering time (or mean free path) is shorter than all other relevant times (or mean free 
paths), the ensemble-averaged particle distribution function can be assumed to be spatially uniform 
and isotropic in p, namely (f(x,p,t)) = (f(p 



Tsvtovich 


1977; 


Melrose 


1980; 


Schlickeiser 


2002) 



The resulting momentum diffusion equation describing evolution of the particle distribution 
can be written as 



d_ 
dt 



1 d 



p 2 D(p)-^(f(p,t)) 



(1) 



where the momentum diffusion coefficients D{p) approximates the rate of interaction with fluctu- 
ating electromagnetic field. Several other terms representing physical process that may influence 
evolution of the particle energy spectrum can be added to the diffusion equation ([1]). In particular, 
one can include continuous energy gains and losses due to direct acceleration (e.g., by shocks) and 
radiative cooling. Furthermore, if the diffusion of particles out of the turbulent region is approxi- 
mated by a catastrophic escape rate (or time t esc ), and if there is a source term Q(p, t) representing 
particle injection into the system, then the spatially integrated (over the turbulent region) one- 
dimensional particle m omentum distribution, n(p,t) = Airp 2 (f(p, i)), is obtained from (see, e.g., 
Petrosian k Liulhooi ) 



dn(p, t) d 
dt dp 



D(p) 



dn(p, t) 
dp 



d_ 

dp 



2 D{p) 
p 



+ (p) n(p, t) 



n(p, t) 



+ Q(p,t). (2) 



Let us further assume presence of an isotropic Alfvenic turbulence described by the one- 
dimensional power spectrum W(k) oc k~ q with 1 < q < 2 in a finite wavevector range k\ < 
k < &2> such that the turbulence energy density f* 2 dkW(k) = (5B) 2 /8 7r is small compared with 
the 'unperturbed' magnetic field energy density, £ = (5 B) 2 / B% < 1. The momentum diffusion 



coefficient in equ ations (PE]) can be then evaluated (e.g., iMelrosd Il968l ; iKulsrud &; Pearce 
Schlickeiser! Il989l ) as 



1969 



D(p) 



— r oc p 



(3) 



2 The Fokker-Planck equation can be also derived straight from the definition of the function f(x,p,t), assuming 
that the interaction of the particles with turbulent waves is a Markov process in which every interaction (collision) 
changes the parti cle energy only by a small amount, and that the recoil of the turbulent modes during the collision 
can be neglected (jBlandford fc Eichlerlll987 ). 
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where A2 = 2n/ki is the maximum wavelength of the Alfven modes, v\ = Pa c is the Alfven 
velocity, and r g = pc/eBo is the gyroradius of ultrarelativistic particles of interest here . Similar 



formulae can be derived for the case of fa st magnetosonic modes (e.g., iKulsrud &: Ferraril 11971 



Achterberg||l98ll ; ISchlickeiser Millerlll998l ) . This allows one to find the characteristic acceleration 
timescale due to stochastic particle-wave interactions, i acc = p 2 /D(p) oc p 2 ~ q / ' (3\. Similarly, the 
escape timescale due to particle diffusion from the system of spatial scale L can be evaluated as 

''esc — 

oc p q , where the spatial diffusion coefficient km = (1/3) c A is given by the appropriate 
particle mean free path, A r* (1/3) C _1 r g (A 2 / r g ) q oc p 9 , that can be fo und from the standard 
relation D(p) = (1/3) (3\ p 2 c/A (for more details see, e.g.. ISchlickeiserl |2002| ) . 



For convenience we define the dimensionless momentum variable x — p/P0i where po is some 
chosen (e.g., injection) particle momentum. With this, the (stochastic) acceleration and escape 
timescales can be written as 



■X 



2 q , where T a 



A 2 



■X 



q 2 , where r P . 



((3{c \eB \ 2 
9L 2 ( ( Poc 



A2 c V ei?o A 



2-q 



q-2 



(4) 



Hereafter we also consider regular energy changes, strictly energy losses, being an arbitrary function 
of the particle energy as given by the appropriate timescale ti oss = hossip)-, namely (p) = —(p/ti oss ). 
We define further r = t/r acc , N(x,t) = p n(p,t)V, and Q(x,t) = r SLCC p Q(p,t)V , where V = 
J d?x is the volume of the system. With such, the momentum diffusion equation ([2]) reads as 



dN 

~B7 



d_ 

dx 



X' 



dN 
dx. 



X$ X )N] -ex 2 - q N + Q. 



or, in its steady-state {dN/dr = 0) form, as 



d_ 

dx 



X' 



dN 
dx 



±.[{2tT1-x* x ) N}-ex 2 - q N 



-Q. 



In the above, we have introduced 



n ''"acc 1 T~i 

v Y = ^— and e 

X iioss(x) 



acc 
''"esc 



(5) 



(6) 



(7) 



Some specific solutions to the equation (0) were presented in the literature. Majority of 
investigations concentrated on the 'hard-sphere approximation' with q = 2, i.e. with the mean 
free path for particle- wave interaction independent of particle energy (A = £ A2/3; 'classical' Fermi- 
II process). It was found, that in the absence of regular energy losses (i? x = 0), the steady-state 
solution of equation © with the source term Q{x) oc 5{x — Xinj), where 5 (x) is the D i rac delta, is 



of a power-law form N(x > Y j n j) oc x " with a = — (1/2) + [(9/4) + e] 1 / 2 (|Davislll956l ; lAchterbere 



19791 ; iPark Sz Petrosianlll9 95). N ote, that for e <C 1 this can be approximated by a ~ l+e/3, which 
the original result obtained by Fermi ( 19491 ). In addition, with the increasing escape timescale, 



is 
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e — ► 0, the steady-state solution approaches N{\ > Xinj) oc x 1 - This agrees with the general 
finding that for the range 1 < q < 2 and the same injection conditions the steady-state particle 
energy distribution implied by the equation ([6]) is N(x > Xinj ) oc X l ~ q : as long as the regular 



energy changes and particle escape can be neglected (i? x = e = 0; lLacombejll979l ; iBorovskv Eilek 



19861 ; iDroge fc Schrickeiserlll98q ; iBecker et al.ll2006l ). The w hole energy range < \ < oo with the 
appropriate (singular) boundary conditions is considered in lPark &: Petrosianl (j!995l ) . 



The analytic investigations of the momentum diffusion equation ([5]) in the q = 2 limit including 
the radiative cooling have concentrated on t he synchrotron-type losses ~d x oc x ( see > however, 
Schlickeiser et al.l 119871 ; ISteinacker et al.lll988l ). The e xtended dis c ussio n on the time-dependent 
evolution for such a case (equation [5]) was presented by iKardashevi ()1962). As for the steady-state 
solution (equation [6]) , it was found that with Q(x) oc 5(x — Xinj) and the range < x < 00 



N(x > Xinj) oc x 



cr+1 



x 

Xeq JJ 



a 



1, 2c7 + 2, 



X 

Xeq 



(8) 



where a is the energy spectral index introduced above, the equilibrium momentum Xeq is defined 
by the i acc = ii oss co ndition (yie l ding ^ Y = x/Xeq) , and U\a,b,z] is a Tric omi confluent hyper- 

For x < Xeq, i.e. 



geometrical function (j J ones! Il97dl ; ISchlickeiserl Il984l ; iPark &: Petrosian 



199 



for the particle momenta low enough to neglect radiative losses, the above distribution function 
has, as expected, a power-law form N(x > Xinj) oc '. For x ^ Xeq and e <C 1, the particle 
energy spectrum approaches N(x > Xinj) oc x 2 exp (— x/Xeq)- That is, as long as particle escape 
is inefficient, a two component stationary energy distribution is formed: a power-law oc x~ 1 at 
low {x < Xeq) energies, and a pile-up bump ('ultrarelativistic Maxwellian distribution') around 
X ~ Xeq- For shorter escape timescale no pile-up form appears, and the resulting particle spectral 
index depends on the ratio e of the escape and the acceleration time scales. 

In the case of q ^ 2 and synchrotro n- type energy losses $ x oc x, the steady-state solution to the 
equation (El) proyided bylMelrosei (ll969n was qu estioned due to unclear boundary conditions applied 



( Tademaru et al. 



1971 



Park &; Petrosianl Il995l ). The special case of q = 1 with particle escape in- 



cluded (and the infinite energy range < x < °°) was considered further bv lBogdan Schlickeiser 
()1985l ). It was found, that with the injection of the Q(x) oc 6(x — Xinj) type, the steady state 
solution of equation ([6 



id 41 



N(x > Xinj) oc x 2 



\ 



u 



Xeq 
Xcsc 



, 2, 



2 Vxc q y 



(9) 



3 The effects of regu lar energy gains were omitted here for clarity. Jones! 1 197p| ); Schlickeiser ( 1984 ) and 



Park fc Petrosianl ([1995) included in their investigations regular energy gains representing very idealized shock accel- 



eration process. Within the anticipated 'hard-sphere' approximation, these gains were assumed to be characterized 
by the appropriate timescale independent on the particle energy, (p) gain OC p . 



4 T he other (non-synchrotron) radiative losses terms included in the analysis presented bv lBogdan fc Schlickeiser 
(1985) were omitted here for clarity. 
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where the critical escape and equilibrium momenta Xesc and Xcq are denned by the conditions 
t esc = t acc and t acc = fjogg, respectively, yielding e = l/xL: and ^x = x/Xeq- Tnis solution implies 
N(x > Xinj) oc const at low particle momenta for which synchrotron energy losses are negligible 
(x ^ Xeq), independent of the particular value of the escape timescale. At higher particle energies, 



an exponential dependence is expected, N(x > Xinj 
that with an increasing escape timescale this approaches 



(X -^-(Xeq/Xcsc) 2 eX p 



"5 (x/XcqY 



Note, 



X exp 



-\ (x/XeqY 



3. Inefficient Particle Escape 

In this section we are interested in steady-state solutions to the momentum diffusion equation 
([6]) in the case of a very inefficient particle escape and a general (i.e., not necessarily synchrotron- 
type) form of the regular energy changes $ x , which is however a continuous function of the particle 
energy. With e = 0, the homogeneous form of this equation can be therefore transformed to the 
self-adjoint form 

G(x)N( X )=0 (10) 



dx 



P(X) ^N( X ) 



with 



P(x) = x q s( x ) 

G(x) = 



d 



2(q-l)x' ! ~ 2 -^(x^ 



s(x) 



s(x) 



X exp 



dx'x n - q # x > 



(11) 



We also restrict the analysis to the finite particle energy range x£ [xi> X2], where < xii X2 < oo- 
The justification for this is that for a finite range of the turbulent wavevectors, say k E [k\, kj\, 
the momentum diffusion coefficient as given in equation ^ is well defined only for a finite range 
of particle energies (momenta). For example, gyroresonant interactions between the particles and 
the Alfenic turbulence require particles' gyroradii comparable to the scale of the interacting modes, 
or fcr g ~ 2tt. Hence, the lower and upper limit of the particle energy range could be chosen as 
Xi = 2-KeBo/kicpQ and X2 = 27reBo/k2cpo, respectivehjf. Since all of the functions P(x)i P'{x)i 
G(x)i S(x) are continuous, and P(x), S(x) are finite and strictly positive in the considered (closed) 
energy interval, the appropriate boundary value problem, 



ai N(xi) + a 2 



dN( X ) 



dx 



0, 



XI 



5 In the case of the magnetosonic-type turbulence, interacting with particles via transit-time damping satisfying 
the Cherenkov condition fcr g < 1, the low energy cut-off in the momentum diffusion coefficient could be chosen to 
be the energy of the particle whose velocity is comparable to the velocity of the fast magnetosonic mode, which is 
~ va for low-/3, or magnetically dominated plasmas. 
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h JV(xa) + b 2 



dN( X ) 



0; 



(12) 



A2 



is regular. If one of these conditions is violated, which is the case for the infinite energy range < 
y < o o, the problem becomes singular, and the extended analysis presented by lPark Petrosian 
(|1995l ) has to be applied. 



The two linearly-independent particular solutions to the homogeneous form of the equation 
IB are 



yi{x) = s 1 (x), 
mix) = s-\ x ) I d x 'x'- g su'), 



(13) 



or any linear combination of these, 



ui(x) = 2/1 (x) + ay 2 (x) , 
u 2 (x) = yi(x) + Py2(x) 



(14) 



(each involving arbitrary multiplicative constants). By imposing the boundary conditions (1121) in 
a form 



Ol «l(Xl) + a 2 
h U2(X2) +h 



dui(x) 



dx 
du 2 (x) 



dx 



Xl 



X 2 



0, 
0. 



(15) 



parameters a and (3 can be determined. With thus constructed particular solutions to the equation 
(fTUl) . one can define the Wronskian iu(x) = u i(x) u ' 2 (x) ~ u i(x) n 2(x)> an d next construct the 
Green's function of the problem, 



G{X,Xmj) 



1 



U\(x) u 2 (Xinj) for Xl < X < Xinj 



xfnj^CXinj) 1 ^l(Xinj)u 2 (x) for Xinj < X < X2 



(16) 



where xi < Xinj < X2- This gives the final solution to the equation 



N( X ) 



\2 



dXinj^(X,Xinj)Q(Xinj)- 



(17) 



Al 



Steady-state solutions exist, however, only for some particular forms of the injection function 
Q{Xi T )- To investigate this issue, and to impose correct boundary conditions for the finite energy 
range Xl < X 5; X2, let us integrate equation ([5j) over the energies and re- write it in a form of the 
continuity equation, 

n\r /-\-' 

(18) 



9r 
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Here M = J* 2 d\N{x) is the total number of particles and the particle flux in the momentum 
space is defined as 

F[N{ X )\ = - x * x ) N - X q ^ ■ (19) 

Note, that with the particular solutions u\(x) and «2(x) given in (fTil) . one has 

H u i(x)\ = -a, 

HMx)} = -P, (20) 

independent of the momentum x ° r of the particular form of the direct energy losses function 

Let us comment in this context on the 'zero-flux' boundary conditions of the type (|12p . namely 
IF\ = IF\ Vn = 0. These, with equations (fT5l) and (|20|) . imply a = ft = 0, i.e., ui(x) = U2(x)- 

XI X2 

In other words, one particular solution yi(x) satisfies the 'no-flux' boundary condition of the ho- 
mogeneous form of the equation (flQ|) for both xi an d X2- I n such a case, the steady-state solu- 
tion can be constructed using the function yi(x) only if it is orthogonal to the source function, 
f^ 2 dx yi(x) Q(x) = 0- This condition, for any n on-zero partic l e iniection and m(x) = S~ 1 (x) as 



given in the equation (jllf) . cannot be fulfilled (cf. IMelrosdll969l ; iTademaru et al.lll97ll ). 'Zero-flux' 



boundary conditions for non- vanishing Q(x) can be instead imposed if the particle injection is 
balanced by the particle escape from the system (see § 3 below). 

In the case of no particle escape, with the stationary injection such that J^ 2 dx Q(x) = A and 
with the direct (radiative) energy losses $ x 7^ 0, the boundary conditions can be chosen as 

-F\ XI =A, and ^| X2 =0, (21) 

which give a = A and (3 = 0, and correspond to the conservation of the total number of particles 
within the energy range [xi>X2]- Let us justify this choice by noting that the radiative losses 
processes, unlike momentum diffusion strictly related to the turbulence spectrum, is well defined for 
particle momenta x < Xi an d X > Xi- Hence, with non- vanishing radiative losses (as expected for 
ultrarelativistic particles considered here), no flux of particles in the momentum space through the 
maximum value X2 toward higher energies is possible (radiative losses in the absence of stochastic 
acceleration will always prevent from presence of particle flux above X2)- F° r the same reason, there 
is always a possibility for a non-zero particle flux toward lower energies through the xi point, since 
the stochastic acceleration timescale, even if being an increasing function of the particle energy, is 
always finite at xi > 0- Note in this context, that the particle flux at xi implied by the chosen 
boundary conditions (I2T1) must be negative, F\ < 0. That is, there is a continuous flux of particles 
through the xi point from high to low energies, which — in the absence of particle catastrophic 
escape from the system — balances particle injection. With these, one can find the Green's function 
as 

S^ix) U~ x + J*dx'x'- q S(x')) for Xi<X<Xinj , 
<?(x,Xinj)| loss = { u A, , A (22) 

S x (x) [A + J d X X q S{x )) for Xinj < X < X2 
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where S(x), introduced in the equation (fTT|) . can be re-written as 

" x d X ' Wx' 



S(x)=X exp 



X! kossiX 1 ) 



(23) 



3.1. Synchrotron Energy Losses 



As an example let us consider synchrotron energy losses of ultrarelativistic electrons, which 
are characterized by the timescale 



tsyn — ''"syn X 



with Tf 



syn 



67r m,gC 2 
a T p B% 



(24) 



(e.g. , iBlumenthal fc GouldlflQToh , and which define the equilibrium momentum %eq — (Tsyn/^acc) 
through the condition i acc = i syn , yelding •& = x/Xeq q ■ The Green's function (f22|) adopts then the 
form 



.3-9 



£(x,Xinj)| syn = X 2 e 3 ~ q 



Xeq, 



A 



+ 



min[Xinj, X] 



dx X 



t \ 3 — q 



(25) 



g 3-9 V Xcq 



q) 



(3 - g)V(3-«) 



1 + g (min[xinj, x]/Xeq) 3 9 (Xi/Xeq) 3 9 
3 — g ' 3 — g 3 — g 



where r[a, 21,22] is generalized incomplete Gamma function. By expressing the above solution in 
terms of Kummer confluent hyper geometrical f unctions M\a,b,z] using the identity T[a, £1,22] = 
a -1 Z2 M[a, 1 + a, — Z2] — a" 1 zf M[a, 1 + a, —z±] (jAbramowitz fc Stegunlll964! ). assuming xi <C Xeq, 
and noting that M[a, b, z] ~ 1 for 2 — ► 0, one can rewrite it further as 



£(x,Xinj)L n ~ X e ~* 



9 - — I 

z " 3-q Vxeq 



,3-9 



syn 

I A 1 + g 



(26) 



min(xinj, X) 



-1-9 



1 + 9 



M 



1 + g 2-2g 1 /min[x in j, x]V 9 



3 — g' 3 — g ' 3 



Xeq 



Finally, noting that M[a,b,z] ~ T(b)e z z a b /T(a) for z — > 00, and neglecting the A 1 term, one 
finds a rough approximation 



£(X,Xinj)| 



syn 



J_ V "1-9 V 2 -^(X/Xcq) 3 9 
l+q Al A e 

Xcq 9 X- 2 
Xeq^Xrnf^^ 7 ^) 3 " 9 



for min(xinj,x) ^ Xcq 

for Xeq < X < Xinj 



X 2 e (X/X ° q)3 9 for Xeq < Xinj < X 



Hence, as long as min[xi n j, x] < Xeq, one has G(x, Xinj)| syn oc x 2 exp --^ (x/Xeq) 



\3-<j 



3-9 



(27) 
If, how- 



ever, min[xinj) x] > Xeq, the Green's function retains the spectral shape oc x 2 exp —3+^ (x/Xeq) 3 9 
for Xinj < X, while is of a power-law form Q(x, Xinj)| syn oc x~ 2 for X < Xinj- 



- 11 - 




Fig. 1. — Upper panel: Stochastic acceleration timescales for fixed plasma parameters (Bq, C,, (3a, 
Xi) X2) but different turbulence energy index: q = 2 (thick solid lines), q = 5/3 (thick dashed 
lines), and q = 1 (thick dotted lines). Thin solid line denotes radiative (synchrotron) energy losses 
timescale considered. Lower panel: Particle spectra resulting from joint stochastic acceleration and 
radiative (synchrotron) energy losses specified in the upper panel. The spectra correspond to the 
monoenergetic injection Q(x) °t <^(x — 1) with fixed J dpQ(p), and no particle escape. Thin solid 
line denotes particle spectrum expected for the same injection and cooling conditions, but with the 
momentum diffusion effects neglected, N(x)- 
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In Figures ([T|) and ([2]) we plot examples of particle spectra obtained from the above solution 
for the system with fixed plasma parameters (Bo, (, /3a, Xi, X2) but different turbulence energy 
indices: q = 2 ('hard-sphere' approximation; thick solid lines in the figures), q = 5/3 (Kolmogorov- 
type turbulence; thick dashed lines), and q = 1 (Bohm limit; thick dotted lines). As for the source 
function, we consider two different forms, namely Q(x) °c <5(x _ 1) i n Figure[T]and Q(x) on <5(%— 10 ) 
in Figure O with the normalizations given in both cases by the same fixed J dpQ(p). The emerging 
particle spectra are compared with the ones expected for the same injection and cooling conditions, 
but with the momentum diffusion neglected, N(x)- Such a steady-state electron distribution can 
be found from the appropriate equation 



d_ 

dx 



X$ X N( X ) +Q(x) = o 



(see equation [5]) , for which one has the straightforward solution (IKardashevI 1 1 962 ) 



N(x) 



X$ 



x Jx 



Q(Xinj)dx 



(28) 



(29) 



(thin solid lines in the lower panels of Figures [TJ{2j) . 

As shown in the figures and follows directly from the obtained solution [2514271 joint stochastic 
acceleration and radiative (synchrotron-type) loss processes, in the absence of particle escape, tend 



to establish N(x) oc x 2 ex P 



(x/Xcq) 



3-q 



spectra independent of the energy of the injected 



3-g VA/ Aeqy 

particles and the form of the source function as long as it has a narrow distribution. Moreover, 
for x ^ Xeq the turbulence energy index q does not influence the spectral shape of the electron 
energy distribution. Instead — with fixed normalization of the injection function Q(p) and fixed 
plasma parameters (including magnetic field intensities Bq and £) — turbulence power-law slope q 
determines (i) the equilibrium momentum Xeq, (h) normalization of the electron energy distribution, 
and (hi) the spectral shape of the particle distribution for x > Xeq- in particular, flatter turbulent 
spectrum leads to higher value of the equilibrium momentum x eq , lower normalization of N(x), and 
steeper exponential cut-off at x > Xeq- Note also, that if particles with momenta xinj 3> Xeq are 
being injected to the system, the resulting electron energy distribution may adopt the 'standard' 
form of the synchrotron-cooled source function (|29p at highest momenta Xeq ^ X < Xmj (e-g-, 
oc x~ 2 f° r the Q(x) oc 5(x — 10 7 ) injection in Figure [2]). 



3.2. Inverse- Compton Energy Losses and the Klein-Nishina Effects 

Let us now investigate the effects of the inverse-Compton (IC) radiative energy losses in the 
presence of a turbulent particle acceleration. At low energies when the Klein-Nishina (KN) effects 
are negligible the IC case is identical to the synchrotron case with the magnetic energy density 
B 2 /8tt replaced by the photon energy density u p h- The two cases differ when KN effects become 
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important at high energies. To include these effects we approximate the radiative loss timescale as 

1.5 

Xcr J 



kc = nc x 



-i 



i + JL 



where tic 



ia T p u ph 



and Xc 



m e c 
4R)£o 



(30) 



Here the radiation field involved in the IC scattering was assumed to be monoenergetic, with the 
total energy density u pn and the dimensionless (i.e., expressed in the electron mass units) photon 
energy en- The above formula properly takes into account KN effect up to energies \ — 10 Xcr 
(jModerski et al.1 120051 ). Clearly, as long as q < 1.5, balance between acceleration and cooling 
timescales takes place at one particular momentum Xeq = m ax(xTh> Xkn)> depending on weather 
energy losses dominate over acceleration in the Thomson regime, Xeq = XTh = (^c/Tacc) 1 ^ 3-9 ^, 
or in the KN regime, Xcq = Xkn = Xtii ^ Xcr 1 ' 5 5 • For q > 1.5, there may be instead 

two equilibrium momenta for a given one acceleration timescale, Xeq, l = XTh and Xeq, 2 = Xkn> 
or no equilibrium momentum at all, if i acc < t lc within the whole considered range x < 10 4 Xcr- 
Finally, for q = 1.5 (that corresponds to the Kraichnan turbulence), the ratio between IC/KN and 



oc X 



9-1 



X 



acceleration timescales is energy-independent, since both t 
kc(x > Xcr) oc X 1/2 - 

Assuming hereafter q ^ 3/2, one can find from the equation f|23|) that 



l / 2 and, as given in ([30 



Six) = X 2 exp 



1 



3- q \XtJ 



3 - q, 4 - q, 



X 

Xcr 



(31) 



where F[a, b, c, z] is Gauss hypergeometric function. This gives the Green's function 



£(X,Xinj)|- c ' " ~ \~ '^P 




miia(xinj,x) 



X 



Q VXTh 



3-g 



dx' X 



l-2-q 



1 



exp 



X 



3-g VXTh 



3-g 



F 



2 



X 

Xcr 



4X 
, - q, 

Xcr 



(32) 



Below we discuss some properties of the obtained solution by expanding the Gauss hypergeometric 
functions as F[a,b,c,z] ~ 1 for z — ► 0, and F[a, b,c, z] ~ [r(c) T(b — a) /Yip) T(c — a)] (—z)~ a + 
[r(c) r(o - b) /Via) Tic- b)] i~z)~ h for z -> oo. 



Let us consider first the case of a low-energy injection, such that Xinj < rniri(XTh, Xc 
Green's function (|32h can be then approximated roughly by 



Six, Xinj) | 



iC|Xinj< 



i x -i-9 e - ^ r(3- g ) r(,-1.5) (Xc r /XT h ) 3 ^ ^2 g - 



(x/Xkn) 



1.5 — q- 



1+9 



The 

for x < Xcr 

for x > Xcr 
(33) 



For x < Xcr the Green's function has the same form as the synchrotron case, which is expected 
in the Thomson regime. However, x > Xer> the KN effects modify the high-energy segment of the 
particle energy distribution. In particular, for q < 1.5 (e.g. q = 1 in Fig. [3]) the spectrum is always 




Fig. 3. — Upper panel: Stochastic acceleration timescales for fixed q = 1 and different plasma 
parameters (thick solid and dashed lines). Thin solid line denotes inverse-Compton energy losses 
timescale considered with the assumed Xcr = 10 4 . Lower panel: Particle spectra resulting from 
joint stochastic acceleration and inverse-Compton energy losses specified in the upper panel. The 
spectra correspond to the monoenergetic injection Q(x) 5{\ — 1) with fixed J dpQ(p), and no 
particle escape. Thin solid line denotes particle spectrum expected for the same injection and 
cooling conditions, but with the momentum diffusion effects neglected, N(x)- 
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log x 

4. — The same as Figure (|3|) except for q = 2. 



-17- 



of a 'single-bump' form, possessing either a sharp or a smooth exponential cut-off depending on 
whether we are in the Thomson or KN cooling regime, respectively. On the other hand, with q > 1.5 
(e.g. q = 2 in Fig. 2]) the acceleration and loss timescales can be equal at two different energies, in 
which case the particle spectra become concave, flattening smoothly from the exponential decrease 



oc x exp 



3^ (x/XTh) 



3-9 



at XTh < X < Xcr to the asymptoticaly approached oc x 2 continuum 
at X > Xkn- Such spectra are shown in Figure (OHJ) for q = 1 and q = 2, respectively, assuming 
monoenergetic injection with Xinj = 1, Xi = 10~ 2 , Xkn = 10 4 , and X2 = 10 8 . In each figure we 
use two different acceleration timescales (thick solid and dashed lines), but the radiative losses 
timescale, ti OS s 5 as well as the normalization of the injection function, J dpQ(p), are kept constant. 
The emerging spectra are compared with the electron energy distribution N(x) corresponding to 
the same injection and cooling conditions, but with the momentum diffusion neglected (equations 
thin solid lines in the lower panels of the figures). 

In the case of q ^ 1.5 and high-energy injection Xinj > Xcr> the appropriate Green's fuction 



retains again familiar shape <5(x,Xinj 



19^1-5 

'ic;xmj> 



— l— g 2 - 
Xi X e 



i 

3-9 



(x/xt) 



3-9 



at low particle mo- 



menta x < Xcr- And again, at % > Xcv significant deviations from such a form may be observed, as 
follows from the approximate form of the Green's function 



^(x,Xinj)|^; 5 

min(xinj,x) 



dx X 



X 2 e" 



t-2-q 



^(X/XKN) 1 ' 5 -^- 



-^r(3- ? )r(<z-i.5) (xcr/xTb) 3 -" x 



(34) 



exp 



XI 



X' 



q vxTh 



3-9 



F 



<1, 4- q, 



X 

Xcr 



for x > Xc 



(see equation 1321 with the A -1 term neglected). The resulting particle spectra are plotted in Figures 
(|5][6]), where we consider two limiting cases of q = 1 and q = 2, and assume monoenergetic injection 
Q(x) °c 5{x ~ Xinj) with Xinj = 10 7 . All the other parameters are fixed as before. As shown, 
in addition to the spectral features discussed in the previous paragraph for the case of a low- 
energy injection (Figures [3JH]) , the radiatively-cooled continuum may be observed at high particles 
energies x < Xinj, depending on the efficiency of the acceleration process. The KN effects manifest 
thereby by means of a characteristic spectral flattening over the 'standard' power-law form oc x -2 > 
obviously only within the momentum range Xcr < X < Xinj) i n agreement with the appropriate 
N(x) distribution (thin solid lines in the lower panels of Figures [M6|). Such a feature, being a 
direct result of a dominant IC /KN-regime radiative cooling with the momentum diffusion effect s 
negligible, was di s cusse d previously by, e.g., iKusunose k, Takaharal (|2005l ); iModerski et al.l (|2005l ); 
Manolakou et al.1 ((2003). 



Finally, for completeness we note that with q = 1.5 one can solve equation (|23p to obtain 
S(x) = X~ 2 exp { 2 ( — ) ( ArcSinhy 7 x/Xc 



Vx/xc 



V 1 + (X/Xcr) 



(35) 



This reduces to S(x) ~ X 2 ex P | (x/XTh) 3 ^ 2 for x < Xcr, and can be approximated by S(x) ~ 
0.54x _1 Xcr 1 f° r X > Xcr- The resulting particle spectra, shown in Figure ([7]) for the case of a 
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log x 



Fig. 5. — Upper panel: Stochastic acceleration timescales for fixed q = 1 and different plasma 
parameters (thick solid and dashed lines). Thin solid line denotes inverse-Compton energy losses 
timescale considered with the assumed Xcr = 10 4 . Lower panel: Particle spectra resulting from 
joint stochastic acceleration and inverse-Compton energy losses specified in the upper panel. The 
spectra correspond to the monoenergetic injection Q(x) 5(x — 10 7 ) with fixed j dpQ(p), and 
no particle escape. Thin solid line denotes particle spectrum expected for the same injection and 
cooling conditions, but with the momentum diffusion effects neglected, N(x)- 
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Fig. 7. — Upper panel: Stochastic acceleration timescales for fixed q = 3/2 and different plasma 
parameters (thick solid, dashed, and dotted lines). Thin solid line denotes inverse-Compton energy 
losses timescale considered with the assumed Xcr = 10 4 . Lower panel: Particle spectra resulting 
from joint stochastic acceleration and inverse-Compton energy losses specified in the upper panel. 
The spectra correspond to the monoenergetic injection Q(x) oc 5(x — 1) with fixed f dpQ(p), and 
no particle escape. Thin solid line denotes particle spectrum expected for the same injection and 
cooling conditions, but with the momentum diffusion effects neglected, N(x)- 
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low-energy injection Q(x) oc $(x ~ 1)> are therefore iV(x < Xcr) oc x 2 exp — | (x/XTh 



^3/2 



at low 



momenta, or of the power-law form N(x > % 
are important. Here a' = t . (x>x ) — ^ = 



1.5 

Acr 



oc X 
2. 



at higher momenta where the KN effects 



3.3. Bremsstrahlung and Coulomb Energy Losses 

At high densities or low magnetic field (in general low Alfven velocities) electron-electron 
and electron-ion interactions become important. These result in an elastic loss due to Coulomb 
collisions or radiative loss via bremsstrahlung. At low energies the bremsstrahlung loss rate is 
negligible when compare d to the Coulomb loss rate, which is independent of energy for relativistic 
charge particles (see e.g. |Petrosianlll973l . l200ll ). However, since the bremsstrahlung rate increases 



nearly linearly with energy, above some critical energy bremsstarahlung becomes dominant. The 
time scales associated with these processes approximately are 



Wi = Tcoul X, where 

Tcoul 



m e c 3cJThC^g In A 



and 



tbr 



Tbr 



where Tb r 



7T 



(36) 



(37) 



3 af s o"Th cre g 

Here n g is the background plasma density, the Coulomb logarithm In A varies from 10 to 40 for 
variety of astrophysical plasma, otf s = 1/137 is the fine structure constant, and the bremsstrahlung 
rate includes electron-ion and electron-electron bramsstrahlung, a nd assumes completely unsc reened 
limit with approximately 10% (fully ionized) helium abundance (IBlumenthal Gouldlll970l ), The 
time scales are equal at energy pcoul = tt In Am e c/(2af s ). At higher energies the bremsstrahlung 
loss becomes unimportant compared to the synchrotron or IC losses. For example, the synchrotron 
loss becomes equal to and exceeds the bremsstrahlung loss at electron momenta p > Pbrem = 
(m e /m p )(a{ s / m e c so that for bremsstrahlung to be at all important we need 1000 < p/ (m e c) < 
10~ 5 /3^ 2 , requiring f3\ < 0.003. Below we investigate in some details stochastic acceleration for 
the conditions when the Coulomb and bremsstrahlung processes are the dominant loss processes. 

At low energies, p < Pcoulj Coulomb collision dominate. If po ^ m e c then in the range 
m e c <C p -C PCoul and for q > 1, the appropriate Green's function becomes 



£(X,Xinj)| coul = X 



9>1 _ ,,2 1-9 



, Xcq 



, Xeq , 



A 



+ 



min[xinj,x] 



dx X 



v Xeq 



Xcq 



XI 



X e 



(1 - qf/^-l) 



1 + 9 



(min[xinj, Xl/Xeq) 1 q (Xl/Xc V ' ' 



1 



1 



(see equations [22T- I23]) . where the equilibrium momentum Xeq = (Tcoui/^acc) 1 ^ 1 ^ is defined by the 
tacc = icoul condition, yielding t? x = Xcq 1 /x- Note that since q > 1 are considered, the acceleration 



(38) 
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timescale is longer than the Coulomb interactions timescale for x < Xeq- Thus, in the case of a 
low-energy particle injection with Xmj < Xeq ; the emerging particle spectra are of the 'cooled' form 
N(x) = -^(x) ^ const (see equation [29] with $ x oc x )■ If> however, higher-energy particles are 
injected to the system, an additional flat-spectrum component N(x) oc x 2 is formed at x > Xeq- 

Let us finally note, that pure Coulomb energy losses and the Bohm limit q = 1 correspond 
to the situation when i? x = const, and hence S(x) = x~ 2+ ^ cc ^ Tcou ^ ■ The Green's function (|22p 
adopts then the form 

„ =1 n rar.n / 1 /-m^IXinj , X] ; , 

^(x,xmj)lcoui = x Tcoul ( A + y x dx x 

1 -V J Xi for T acc /r cou i < 2 

a' ^ mm' (xinj, Xj for Tacc/Tcoul > 2 

where a' = — 2. Hence, if only r acc < 2r cou i, a power-law particle energy distribution N(y) oc 

Tcoul 

X~ a forms, with — 2 < a 1 < 0. For any longer acceleration timescale, r acc > 2t cou i, and for 
the source function Q(x) oc 5(x — Xinj)> the emerging electron spectra are N(x) oc const for 
X < Yinii and iV(x) a Y~ g ' w ith cr' > for x > Xinj- This is consistent with the solution found by 



Bogdan Schlickeiserl ([1985), who considered also synchrotron emission and finite escape timescale 
in addition to the Coulomb energy losses of ultrarelativistic electrons interacting with flat-spectrum 
turbulence q = 1. 

At higher energies and in the range Pcoul ^ P Pbrem bremsstrahlung loss is the dominant 
process and the equilibrium momentum defined by the condition t acc = ibrem f° r Q < 2 becomes 
Xeq = (Tbrem/Vacc) 1 ^ 2- ^, yielding d x = Xeq ?) • Hence, the Green's function ([22]) is 

n 1 ( * \ 2 ~ q ( 1 /•min[ X inj,x] i ( x ' \ 2 ~ q \ 

^(x,Xinj)|^ e 2 m = X 2 e-^V^J +/ dx'x'-^e^^) U (40) 




~ 2 -^f^) 2 ' 9 Xe-q 1 -^-!)^ 2 - 9 ^ ^ 1 + g (min[ Xinj , X}/X cq ) 2 - Q (Xl/Xcq^ 9 
~ X 6 ( 2 -g)3/(2-g) 

(equations [221— I23j) . In other words, for any injection conditions the expected electron energy 



distribution is of the N(x) oc x 2 ex P 



\2-<? 



form, except for the case when high energy 



~ q (x/Xcq, 

particles with Xinj > Xeq are injected to the system. "Such high energy particles subjected to the 
bremsstrahlung energy losses form then an additional 'cooled' high-energy power-law tail N{x) oc 
X _1 in the momentum range between Xeq and Xinj) in agreement with the appropriate form of iV(x) 
with -& x = const (see equation [29]) . 

The situation changes for q = 2, since both the acceleration and cooling timescales are now 
independent of electrons' energy. In this case S(x) = x _2+ ^ Tacc ^ Tbrcm \ and the Green's function (|22p 
adopts the form 

^(X,Xinj)lbrem = X brcm T + / dX X ~ 
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1 x -a' f Xl lW for Tacc/r b rcm < 3 

l-o-' 1 min" 1+,T '(xmj, x) for T acc /r brem > 3 



where cr' = r ^ ac ^ —2. This is consistent with the appropriate Green's function found by lSchlickeiser et al 



(|1987l ) who, in a framework of the 'hard-sphere' approximation q = 2, considered also synchrotron 
emission and particle escape in addition to the bremsstrahlung radiation. The solution (|4ip implies 
that within the whole energy range the expected electron energy distribution is of the power-law 
form N(x) °c X with the power-law index —2 < a 1 < 1. For any longer acceleration timescale, 
Tacc > 3Tb re m> and monoenergetic injection Q(x) oc 5(\ — Xinj), the emerging electron spectra are 
expected to be of the N(x) oc x~ l form for x < Xmj; while N(x) oc with a' > 1 for x > Xiny 



4. Efficient Particle Escape 



In this section we investigate steady-state solutions to the momentum diffusion equation of 
radiating ultrarelativistic particles with a finite escape timescale (equation [6]) . Our analytical ap- 
proach force us to consider only the limiting cases of turbulent spectral indices q = 2 or q = 1, as well 
as to restrict the analysis of radiative losses to the synchrotron and/or IC-Thompson regime pro- 
cesses, (t? x oc x)- We note that the global approximation to the solution of the momentum diffusion 
equation not necessarily restricted to some particular values of the q parameter, with the regular 



energy losses and particle escape terms included, were studied by iGallegos-Cruz Perez-Peraza 



(|1995l ) by using the WKBJ method. Just us before, we consider finite energy range of particles 
undergoing momentum diffusion, < Xi> X2 < 00 ; strictly related to the finite wavelength range 
of interacting turbulent modes. We construct the Green's function accordingly to the procedure 
outlined in the previous section § 3, with addition of the escape term (e ^ 0) and with a different 
boundary conditions. Specifically, we change equation ([5]) to 



cW 
dr 



XI 



X2 



XI 



X-2 



d X Q(x,r)-e I d XX 2 ~ q N(x) 



(42) 



A'l 



where the particle flux in the momentum space J-[N{x)] is defined in the same way as previously 
(equation [T9j) . As evident, the no-flux boundary conditions, .F[./V(xi)] = J r [N(x2)] = 0, and 
conservation of total number of particles, dM/dr = 0, (within the energy range [xi>X2]) implies 
that the particle injection is completely balanced by the particle escape. We will assume this 
to be the case in this section. Physical realization of these would imply presence of an another 
efficient yet unspecified acceleration process operating at x < Xl> which prevent negative particle 
momentum flux through the xi boundary. As shown below, the solutions we obtain agree with the 
ones discussed in the literature for singular boundary cond i tions for the infinite m omentum range 
(|Joneslll970l ; ISchlickeiserl Il984l ; iBogdan &; Schlickeiserlll985l ; iPark &; Petrosianlll995l ). as long as we 
are dealing with particle momenta % ^ Xi anci X ^ X2- 
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4.1. 'Hard-Sphere' Approximation 



'Hard-Sphere' approximation for the momentum diffusion of ultrarelativistic electrons under- 
going synchrotron energy losses corresponds to the fixed q = 2 and r d x = x/Xcq (see equations H] 
and [24]) . With these, the equation ([6]) adopts the form 



X 2 N"(x) + Xcq 1 X 2 N'(x) + (2 Xcq 1 X - 2 - e) N( X ) = -Q( X ) • 



(43) 



The two linearly-independent particular solutions to the homogeneous form of the above equation 
are 



yiix) 

2/2 (x) 



X 



e Xc< i U 



a- 1, 2cr + 2, — 

Xeq. 



X CT+1 e xe q m 



a -1,2a + 2, 



X 

Xeq 



(44) 



where U[a,b,z] and M[a,b,z] are Tricomi and Kummer confluent hypergeometrical functions, re- 
spectively, and a = —(1/2) + [(9/4) -l-e] 1 / 2 . Introducing next their linear combinations, tii(x) = 
2/1 (x) + a 2/2(x) an d U2(x) = 2/1 (x) + PU2(x)i one may find that the no- flux boundary conditions 
F[u\{xi)\ = J 7 [u2(X2)] = are fulfilled for 



a = (2 + a) 



U[a, 2a + 2, Xl/Xeq] 



M[a, 2a + 2, Xl / Xeq ] 
This gives the Green's function of the problem as 



and /? = (2 + a) 



U[a, 2a + 2, x 2 /Xeq] 
M[a, 2a + 2, X2 /Xcq] 



(45) 



S<X, X^llZ 2 = (« - (i)- 1 X- 2 Xeq 2 ' 7 - 1 x 

x f [2/1 (x) +ay 2 (x)] [2/i(Xinj) + /?2/2(Xinj)] for Xi < X < Xinj 
1 [2/i(Xinj) + ay 2 (Xinj)] [2/1 (x) + 0y 2 {x)} for Xinj < X < X2 



(46) 



In order to investigate the above solution, let us consider first the case X i Xinj <C Xeq <C X 2j 



and use the standard expansion of the confluent hypergeometrical functions: U[a, b, z] 



~ z 



and 



T(b)e z z a - b /T(a) for z -» 00, while C7[a,6,z] ~ T(6 - l);* 1 " 6 /!^) and M[a,b,z] ~ 1 



M[a,6,z] 

for 2 — ► (jAbramowitz &: Stegunlll964l ). In this limit one gets 



1 9=2 

'esc, Xinj < 



~ < 



1 



2ct+1 ^inj 



xr 



1 



2<T+1 
r(cr-l) cr-1 
r(2cr+2) ^inj Xcq 



Xinj X 



~ 2 X 2 e 



-X/Xe- 



for Xi < X < Xinj 

for Xinj < X < Xcq 
for Xcq < X < X2 



(47) 



Thus, by moving the critical momenta xi an d X2 toward and 00, respectively, the resultant Green's 
function approaches asymptotically — as expected — the correspond ing Green's function fo r singu - 
lar boundary conditions obtained by [Jones] (|197Q ) : ISchlickeiserl (|1984l ) and lPark &: Petrosianl (| 19951 ). 
In particular, one can find that with the monoenergetic injection Q(x) °c <5(x ~~ Xinj); the resulting 
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Fig. 8. — 'Hard-sphere approximation' (q = 2): particle spectra resulting from joint stochastic 
acceleration, particle escape, and synchrotron energy losses. The spectra correspond to the mo- 
noenergetic injection Q(x) °c $(x ~ Xinj) with fixed normalization, fixed acceleration and cooling 
rates, but different escape timescales (parameter e = 3, 0.1, 10~ 4 ; dotted, dashed, and solid lines, 
respectively). For illustration, xi = 1CT 2 , Xinj = 1, Xcq = 10 6 , and \2 = 10 8 have been selected. 
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electron energy distribution is then of the form N(x < Xinj) X CT an d N( X > Xinj) oc X _<J U P 
to maximum momentum Xeq- Moreover, for the increasing escape timescale e —* 0, one has cr ~ 1 
and the pile-up bump iV(x) oc x 2 ex P [ — X/Xeq] emerging around x ~ Xeq energies. This is shown 
in Figure ([8]), where we fixed normalization of the monoenergetic injection J dpQ(p), acceleration 
and losses timescales, but varied the escape timescale (e = 3, 0.1, 10~ 4 ; dotted, dashed, and solid 
lines, respectively). For illustration we have selected xi = 1CP 2 , Xinj = 1> Xeq = 10 6 , and X2 = 10 8 . 

When high energy particles are injected to the system, such that xi ^ Xeq *C Xinj *C X2, one 
may find useful asymptotic expansion of the Green's function 



£(X,Xinj)| 



esc, Xinj > 



r(2cr+2) A eq A 



+1 p-x/xc 



X 2 Xeq 



Xi^ Xeq e Xin j /x 



cq ^2 g- 



-X/Xe 



for Xl < X ^ Xeq 
for Xeq < X < Xinj 
for Xinj < X < X2 



(48) 



Xinj) with Xinj > Xeq the resulting electron 
-X/Xeq] for x ^ Xeq- However, within the 



That is, for the monoenergetic injection Q(x) oc 5(x - 
energy distribution is of the form iV(x) oc x a+1 ex P [ 
energy range Xeq < X < Xinj the power-law tail N(x) oc x -2 emerges, representing radiatively 
oc x) cooled high-energy particles injected to the system, undergoing negligible (when compared 
to the energy loss rate) momentum diffusion. At even higher energies, x > Xinj) the particle 
spectrum cuts-off rapidly. This is shown in Figure ([9]), where, as before, we fixed normalization 
of the monoenergetic injection J dpQ(p), acceleration and losses timescales, but varied the escape 
timescale (s = 3, 0.1, 10~ 4 ; dotted, dashed, and solid lines, respectively). For illustration we have 
selected xi = 10~ 2 , Xeq = 10 2 , Xinj = 10 6 , and X2 = 10 8 . Note, that the esape timescale, and hence 
parameter e, influences now the slope and normalization of particle energy distribution only in the 
'low-energy' regime x < Xeq> such that the spectrum approaches oc x 2 for e — > 0. 



4.2. Bohm Limit 



Bohm limit for the momentum diffusion of ultrarelativistic electrons undergoing synchrotron 
energy losses corresponds to q = 1 and $ x = x/x 2 q ( see equations 141 and 124"]) . The difference with 
the 'hard-sphere' approximation is that the balance between acceleration and escape timescales, 
ta.ee = tesc, define now yet another critical energy, Xesc = £ -1 / 2 and equation ([6]) takes the form 

XA r/ '(x) + (Xeq 2 X 2 -l) A r/ (x) + (2Xeq 2 X-Xe" S eX) N( X ) = -Q(x) ■ (49) 

The two linearly-independent particular solutions to the homogeneous form of the above equation 
are 



yi(x) = x e 2 Vxcqy u 
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where rj = ^ (Xeq/Xcsc) 2 - Defining ui(x) = 2/1 (x) + ay 2 (x) and u 2 (x) = 2/1 (x) +PV2(x), one rinds 



where n = gUWXesc) ■ Uennmg ui(x) = 2/1 (X) + "2MXJ and n 2(X) = 2/HXJ + 
that the no-flux boundary conditions .F[ui(xi)] = T\u 2 (X2)] = corresponds to 



a = 2- 
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+ 1, 3, \ (Xl/Xeq) 
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7?+ 1, 3, 5 (Xl/Xeq)' 



and /3 = 2 
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+ 1, 3, \ (X2/Xcq) 
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' + 1, 3, I (X2/Xeq) 



(51) 



This gives the Green's function of the problem as 



x I [2/1 (x) + oty 2 (x)] [2/i(Xinj) + /32/2(Xinj)] for Xi < X < Xinj 
\ bi(Xmj) + a2/2(Xinj)] [2/1 (x) + /32/2(x)] for Xinj < X < X2 



(52) 



Let us consider first the case xi ^ Xinj *C Xcq X2 for which he Green's function of equation 
521) can be then approximated as 



^(X, Xinj) I 
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for Xi < X < Xinj 

for Xinj < X < Xcq 
for Xeq 



(53) 



< X < X2 



Note that, as expected, in the limits xi - ► and X2 - ► 00, the Green's fu nction (1521) approaches 
asymp totically the solution obtained for singular boundary conditions by iBogdan Schlickeiser 
(jl985l ). As shown in Figure ([10]) . for a monoenergetic injection Q(x) oc <5(x — Xinj), the resulting 



electron energy distribution is N(x < Xinj) oc x arid iV(x > Xinj) oc const up to maximum mo- 
mentum Xeq, with the spectral indexes independent of the value of the escape timescale. However, 
for energies near and above Xeq the spectra depend on the value of rj. For rj — > 0, i.e. when the 
escape timescale is large, the familiar bump N(x) oc x 2 exp [— ^(x/Xeq) 2 ] emerges around x ~ Xeq 
energies (solid line). In the opposite case, when 77 > 1 (or Xeq > Xesc), no pile-up bump is present, 
and the electron spectrum cut-offs exponentially at Xesc momenta (dashed and dotted lines). Here, 
as before, we fixed the normalization of the monoenergetic injection J Q(p) dp, and the acceleration 
and loss timescales, but varied the escape timescale such that Xesc = 10 5 , 10 6 , and 10 7 (dotted, 
dashed, and solid lines, respectively). We choose xi = 1CT 2 , Xinj = 1, Xeq = 10 6 , and X2 = 10 8 . 



In the case when xi *C Xeq *C Xinj *C X2 the asymptotic expansion of the Green's function 
52} yields 
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for Xcq < X < Xinj 
Xinj 



(54) 



for Yi n i < X < X2 



Again as above, the spectrum is different in the case of high energy injection. For example, as 
shown in Figure (fTT|) . for the monoenergetic injection Q(x) oc <5(x — Xinj) with Xinj > Xeq the 
resulting electron energy distribution is of the form N(x) oc x 2 ex P [ — ^(x/Xeq) 2 ] for x ^ Xcq, 
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Fig. 10. — Bohm Limit (q = 1): particle spectra resulting from joint stochastic acceleration, 
particle escape, and synchrotron energy losses. The spectra correspond to the monoenergetic 
injection Q(x) °c $(x ~ Xinj) with fixed normalization, fixed acceleration and cooling rates, but 
different escape timescales (critical momenta Xcsc = 10 5 , 10 6 , 10 7 ; dotted, dashed, and solid lines, 
respectively). For illustration, \i = 10 -2 , Xinj = 1, Xeq = 10 6 , and \2 = 10 8 have been selected. 
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while N(x) oc x 2v ~ 2 f° r Xcq <C x < Xinj- It is interesting to note that the Bohm limit case behaves 
differently from the q = 2 case and analogous injection condition. The escape timescale affecs now 
(via the parameter Xeac, ° r v) the normalization of the low-energy {x < Xeq) segment of the particle 
spectrum but not its power-law slope. It determines, on the other hand, the 'radiatively-cooled' 
part of the particle distribution in the range Xcq < X < Xinj; which is, however, very close to the 
standard oc x~ 2 f° r an Y Xcsc 3> Xcq (or w <C 1). Here, as before, we fixed normalization of the 
monoenergetic injection J dpQ(p), and the acceleration and loss timescales, but varied the escape 
timescale such that Xesc = 10 5 , 10 6 , and 10 (dotted, dashed, and solid lines, respectively). Also 
we set xi = !0" 2 , Xinj = 1, Xeq = 10 6 , and X2 = 10 8 . 



5. Emission Spectra 

In the previous sections § 3 and § 4, we showed that stochastic interactions of radiating ultrarel- 
ativistic electrons (Lorentz factors 7 = p/m e c 3> 1) with turbulence characterized by a power-law 
spectrum W(k) oc k~ q result in formation of a 'universal' high-energy electron energy distribution 



n e (-f) = n 7 exp 



I (jlY 



(55) 



as long as particle escape from the system is inefficient and the radiative cooling rate scales with 
some power of electron energy. Here the equilibrium energy 7 eq is defined by the balance between 
acceleration and the energy losses timescales, while the parameter a depends on the dominant 
radiative cooling process and the turbulence spectrum. In particular, for either synchrotron or 
IC/Thomson-regime cooling one has a = 3 — q. In the case of dominant IC/KN-regime energy 
losses (with q < 1.5) one has instead a = 1.5 — q. Below we investigate in more details emission 
spectra resulting from such an electron distribution. 



5.1. Synchrotron Emission 

Assuming isotropic distribution of momenta of radiating electrons with energy spectrum n e (7), 
the synchrotron emissivity can be found as 



3v, syn(^J 
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where v c = 3eB / 4irm e c, 
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(57) 



and (z) is a modified Bessel function of the second order (jCrusius fc Schlickeiserlll986l ). Rel- 
atively complicated function (|57p can be instead conveniently approximated by lZ(x) ~ 1.81 x 
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(l.33 + x~ 2 / 3 ) e~ x ( Zirakashvili Aharonianll2007l ). allowing for some analytical investigation 
of the integral (156j) . In particular, one may find that the synchrotron emissivity in a frequency 
range v < u 8yn = ^ c 7^j is of the form j UySy n{v < ^syn) °c i/ 1 / 3 , as expected in the case of a very hard 
(inverted) electron energy distribution at low energies, n e {^ < 7 eq ) oc 7 2 . At higher frequencies, 
however, the synchrotron spectrum steepens. In order to evaluate such a high-frequency spectral 
component, we use the introduced approximation for lZ(x), electron spectrum as given in (|55p . and 
with these we re-write synchrotron emissivity (|56[) as 



ju, syn( zy ) 



1.81^/3e 3 J B 7 3 q ra 



dyg(uj,y) exp[-uh(u,y)] , 



(58) 



47r m e c 2 

where to = u/u syI11 y = 7/7^, g(u, y) = y 2 (1.33 + or 2 / 3 ^ 3 ) -1 / 3 , and h(cu, y) = y' 2 + y a /(u a). 
With such a form it can be noted that for large uj, i.e. for v > v^ yn , the integr al of interest can 



be perform approximately using the steepest descent method (see iPetrosian 

jv,syn(y ^ ^syn) — 
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(59) 



where y± = (2uj) 1 ^ 2+a ^ is a global maximum of h(oj,y), and h"(uj,y) = d 2 h(uj,y)/dy 2 . Thus, the 
high-energy synchrotron component drops much less rapidly than suggested by the emissivity of a 
single electron, lZ{x) oc e~ x . For example, assuming synchrotron (and/or IC/Thompson-regime) 
dominance a = 3 — q, the synchrotron emissivity reads very roughly as 



ju, syn {v > Vi 



syn. 



oc v 1 ! 2 exp 



-1.4 {v/v, 



syn; 
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(61) 



In the case of the IC/KN-regime dominance, a = 1.5 — q, the emerging high-energy exponential cut- 
off in the synchrotron continuum can be even smoother than this, for example j v>S yn {y > v t 



syn. 



OC 



v exp 



-2.9 (u/v, 



s0.2 



syn./ for 9 = 1. These spectra are shown in Figure (|12j) for fixed parameters 
.B, no, and 7 eq , where both integration of the exact form of TZ(x) as given in equation (I57h was 
performed (solid lines), and also approximate formulae following from (|59|) were evaluated for 
comparison (dashed lines). Different cases for the parameter a are considered in the plot, namely 
(a) a = 3 — q with q = 1, (b) a = 3 — q with q = 2, and (c) a = 1.5 — q with g = 1. As shown, 
synchrotron spectra are curved and extend far beyond equilibrium frequency ^ sy n- In the case of 



6 As shown by 
electrons. 



Petrosian 



(|198lr ). the following spectral form is also true for synchrotron emission by semirelativistic 
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Fig. 12. — Synchrotron spectra resulting from the electron energy distribution (I55p for fixed pa- 
rameters B, no, and 7 eq . Solid lines correspond to integration of the exact form of lZ(x) as given in 
equation (|57p . while dashed lines to the approximate formulae following from (|59p . Different cases 
for the parameter a are considered in the plot, namely (a) a = 3 — q with q = 1, (b) a = 3 — q with 
q = 2, and (c) a = 1.5 — g with q = 1. 
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the dominant IC/KN-regime cooling with q = 1, the vj v (v) — v synchrotron spectrum peaks around 
~ 10 3 f syn . We emphasize that the approximation (|59p . although obviously not accurate in a range 
v < Vsfsxi works relatively well at higher frequencies, where the standard ^-approximation for the 
synchrotron emissivity, uj VjByn (u) cx [7 s ^(7)] 1/2 , fails. 



5.2. Inverse- Compton Emission 



Let us consider inverse-Compton emission of isotropic electrons up-scattering monoenergetic 
and isotropic photon field with energy density u p h and dimensionless photon energy €q = hvQ/m e c 2 . 
The appropriate emissivity can be then found from 
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167T rripC- 
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ie (l+Vl+Ceeo)- 1 ) T 



^72^^( e ' e 0'7) n e(7) , 



(62) 



where e = hv/m e c 2 , and J7"(e, £0,7) is the IC kernel 



J(e,e ,7) =2JlnJ + J+l-2T 2 + 
(e.g., iBlumenthal Gould! Il97cl ). 
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(63) 



Let us discuss first the case when the KN effects are negligible. The IC kernel can then be 
approximated by J(e,e ,-f) ~ |(l-cj/y 2 ), withy = 7/7 eq , w = e/e ic/Th , and e ic/Th = 4e 7l q which 
is the characteristic energy of soft photon inverse-Compton up-scattered in a Thomson regime by 
electrons with Lorentz factor 7 eq . Hence, with the electron energy distribution of the form (|55p . 
one can find that 



ej e ,ic/Th(e) = ^ co"Th w p h n 7cq dyu? ^1 
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where T\a, z] is incomplete Gamma function. With the expansion T[a, z] ~ T[a] for z —* 
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(jAbramowitz fc Stegunlll964l ). one can approximate further 
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In other words, the IC emissivity at low photon energies is of the form j e ,ic/Th( e < e ic/Th) K e - 
This is the flattest IC/Thomson-regime spectrum, being analogous to the flattest synchrotron one 
ju,syn(v < ^syn) oc u 1 ^ . At higher photon energies, noting that T[a, z] ~ z a ~ l e~ z for z — ► 00, one 
may find instead 
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Fig. 13. — Inverse- Compton spectra produced in the Thomson regime, resulting from the electron 
energy distribution ([55]) for fixed parameters B, uq, and 7 eq . Solid lines correspond to the formulae 
(|64p . and dashed lines to the rough approximation (|66p . Two different parameters a = 3 — q are 
considered in the plot, corresponding to the turbulence energy index q = 1 and q = 2 (cases (a) 
and (b), respectively). 
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Fig. 14. — Inverse- Compton spectra produced in the KN regime, resulting from the electron energy 
distribution ([55]) for fixed parameters B, no, and 7 eq . Solid lines correspond to the exact evaluation 
of the integral (|62p . and dashed lines to the rough approximation (|67p . Different cases for the 
parameter a are considered in the plot, namely (a) a = 3 — q with q = 1, (b) a = 3 — q with q = 2, 
and (c) a = 1.5 — q with q = 1. For illustration 7 cr /7cq = 0.01 has been selected. 
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Therefore, the exponential cut-off of the IC/Thomson-regime component is now steeper than the 
exponential cut-off of the synchrotron component originating from the same particle distribution. 
In particular, with a = 3 — q one gets i e ,ic/Th( e > e ic/Th) <x e ^ 2 ex P[ — \ ( e / e ic/Th)] f° r <7 = 1> while 
Je,ic/Th( e > e k/Th) <x e ex P[ — (eAic/Th) 1 ^ 2 ] for g = 2 (that can be compared with the corresponding 
synchrotron emissivities provided above). These spectra are shown in Figure (|13p for fixed pa- 
rameters B, no, and 7 cq . Here the solid lines correspond to the formulae (|64p . and dashed lines 
to the rough approximation (166p . Two different parameters a = 3 — q are considered in the plot, 
corresponding to the turbulence energy index q = 1 and q = 2 (cases (a) and (b), respectively). 

Finally, we comment on the emission spectra produced in a deep KN regime of the IC scattering, 
i.e. when 7 > 7 cr = l/4eo, by the highest-energy electrons 7 > 7 eq . In such a case, the emissivity 
has to be evaluated by performing the integral (|62p with the exact IC kernel as given in equation 
(|63p . A rather crude approximation for such can be obtained by utilizing the ^-approximation for 
the resulting IC/KN-regime photon energy, namely 6 = 7. In particular, with the electron energy 
distribution as given in (|55|) . and with all the previous assumptions regarding monoenergetic and 
isotropic soft photon field, one finds 



(-)!• < 67 > 

V7cq/ . 

where tic (7) is the inverse-Compton cooling timescale as introduced previously in equation (|30p . 
As shown in Figure [T4"l as a result the IC/KN-regime spectra cut-off sharply above e = 7 eq photon 
energies, imitating exponential cut-off in the energy distribution of radiating particles. Here the 
exact calculations are plotted as solid lines, and rough approximation ([67]) as dashed ones. We fix 
parameters B, no, 7 eq , 7 cr /7eq = 0.01, and again, different cases for the parameter a are considered; 
(a) a = 3 — q with q = 1, (b) a = 3 — q with q = 2, and (c) a = 1.5 — q with q = 1. We also choose 
for illustration 7 cr /7eq = 0.01. 

6. Discussion and Conclusions 

In this paper we study steady-state spectra of ultrarelativistic electrons undergoing momen- 
tum diffusion due to resonant interactions with turbulent MHD waves. We assume a given power 
spectrum W{k) oc k~ q for magnetic turbulence within some finite range of turbulent wavevectors 
k, and consider variety of turbulence spectral indices 1 < q < 2. For example, q = 1 corresponds 
to the 'Bohm limit' of the stochastic acceleration processes, q = 2 represents the 'hard-sphere ap- 
proximation', while q = 5/3 and q = 3/2 to the Kolmogorov or Kreichnan turbulence, respectively. 
Within the anticipated quasilinear approximation for particle-wave interactions, such a turbulent 
spectrum gives the momentum and pitch angle diffusion rates oc p q ~ 2 , or the acceleration and es- 
cape timescales i acc oc p 2 ~ q and t esc oc p q ~ 2 . In the analysis, we also include radiative energy losses, 
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being an arbitrary function of the electrons' energy. In most of the cases, however, or at least in 
some particular energy ranges, the appropriate timescale for the radiative cooling scales simply with 
some power of the particle momentum, ti oss oc p r . For example, r = — 1 corresponds to synchrotron 
or inverse-Compton/Thomson-regime energy losses, r = (roughly) to the bremsstrahlung emis- 
sion, r = +1 (roughly) to the Coulomb interactions of ultrarelativistic electrons, while r = 1/2 may 
conveniently approximate inverse-Compton cooling in the Klein-Nishina regime on monoenergetic 
background soft photon field. 

We find that when the particles are confined to the turbulent acceleration region (t esc — ► oo), 
the resulting steady-state particle spectra (for a finite momentum range of interacting electrons) are 
in general of the modified ultrarelativistic Maxwellian type, n e (p) oc p 2 exp [— - (p/p eq ) a ^, where 
a = 2 — q — r ^ 0. Here p eq is the momentum at which the acceleration and radiative loss timescales 
are equal, i a cc(Peq) = ^ioss(Peq)- This form is independent of the initial energy distribution of the 
electrons as long as this distribution is not very broad and the bulk of initial particles have p < p eq . 
However, if high energy particles with p > p eq are injected to the system, there will be significant 
deviations from this simple form. For example, for a 5-function initial distribution the spectrum 
will have a power-law tail oc p r_1 in addition to the modified Maxwellian bump. Also, if the ratio 
of acceleration and energy losses timescales is independent of the electron energy, in other words, if 
2 — q = r, then the resulting particle spectra are of the form n e (p) oc p~ a ' , where a' = (i a ccAk>ss) — 2. 
Finally, if the particle escape from the acceleration site is finite but still inefficient, a power-law tail 
oc p l ~ q may be present in the momentum range pj n j < p <C Pcqj again in addition to the modified 
Maxwellian component. When the radiative losses timescale is not a simple power-law function of 
the electron energy, the emerging spectra may be of a more complex (e.g., concave) form. 

We also analyze in more details synchrotron and inverse-Compton emission spectra of the 
electrons characterized by the modified ultrarelativistic Maxwellian energy distribution. In order 
to summarize briefly our findings, let us define the critical synchrotron frequency of the electrons 
with the equilibrium Lorentz factor 7 eq = p eq /m e c, namely v syn = (2>eB / A-Km e c) j 2 q , and the critical 
dimensionless energy of the monochromatic (huo = eo m e c 2 ) soft photon field inverse-Compton up- 
scattered (in the Thomson regime) by the 7 eq electrons, e iC/ /xh = ^ e o7eq- With these, one can 
note that the low-frequency synchrotron emissivity is of the form j u ,syn{ v < ^syn) °^ z/ 1 / 3 , as 
expected in the case of a very flat (or inverted) electron energy distribution n e (7 < 7 eq ) oc 7 2 . 
Such flat electron spectra see m to be required to expla in several emission properties of relativistic 
jets in active galactic nuclei (jTsang fc Kirkl l2007al lbl). At higher frequencies, we find a rough 



approximation j v ,syu{ v > u syn) oc i^ 6 a )/( 4+2a ) e xp 
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Thus, the high-energy 



synchrotron component drops much less rapidly than suggested by the emissivity of a single electron 
and the emerging high-frequency tail of the synchrotron spectrum is of a smoothly curved shape 
It is therefore very interesting to note that almost exac tly this kind of curvature is observed at 



synchrotron X-ray frequencie s i n several BL Lac o bjects (jMassaro et al.ll2004l . 120061 : iPerlman et al 



20051 ; iTramacere et al. 
energies. 



2007al fbLI Giebels et al.ll2007l ) , in particular those detected also at TeV photon 
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As for the inverse-Compton emission of ultrarelativistic electrons characterized by the modified 
Maxwellian energy distribution, we find that in the Thomson regime it is of the form J €j i c /Th( e < 

eic/Th) oc e, and j eiic / T h(e > ^ic/Th) oc e (3 ~ a)/2 exp -\ (e/e ic/Th ) a/2 . Both very flat low-energy 
part of this component and also its curved high-energy segment may contribute to the observed 
7-ray emission of some TeV blazars (jKatarzviiski et al.ll2006al : biebels et ap2007l M We also note, 



that the curvature of the high frequency segments of the synchrotron and inverse-Compton spectra, 
event though being produced by the same energy electrons and in the Thomson regime, are different. 
Such a difference is even more pronounce when the Klein-Nishina effects play a role, since in such 
a case an exponential decrease of the high-energy photon spectra is the strongest, j £i i c /KN( e > 
7eq) oc e 7 / 2 exp \—\ (e/7 eq ) a ], imitating exponential cut-off in the energy distribution of radiating 
particles. 
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